On the origin of biquadratic exchange in spin 1 chains 
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One dimensional spin 1 systems may have a rich phase diagram including Haldane gap and dimerized 
phases if the usually very small biquadratic exchange becomes significant. We show that this unlikely 
condition may be fulfilled in electron systems with quasi-degenerate orbitals. This mechanism 
may have been experimentally realized in the spin 1 chain LiVGe2C>6. The implications for the 
exploration of the physics and quantum chemistry of spin 1 chains are discussed. 
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There has been a lot of activity recently on the physics 
of one-dimensional (ID) spin 1 systems, especially after 
the prediction by Haldane that the ID Heisenberg model 
has a spin gap for integer spins [SJ. That prediction has 
been confirmed since then by the observation of a gap 
in many spin 1 chains 0. There has also been a lot of 
progress in the study of the most general Hamiltonian de- 
scribing an isotropic coupling between neighboring spins 
1, namely |3j 

h = y^Hjj, 

(ij) 

Hij — J\Si ■ Sj + J-i(Si ■ Sj) 2 , (1) 

where the sum (ij) is over nearest-neighbor pairs. The 
phase diagram of this model is extremely rich It 
is most easily described using the parametrization J\ = 
Jcos8, J2 = J sin 9. As illustrated in Fig. 1, the system 
has two gapped phases: The Haldane phase for — it/ 4 < 
9 < 7r/4 and a dimerized phase for — 37r/4 < 9 < —ir/4. 
The two phases are connected by a critical point at 9 = 
— 7r/4 for which the model is exactly solvable and the 
spectrum is gapless. In addition, there is a Lifshitz point 
and a disordered point in the Haldane phase, as well as 
a Valence Bond Solid point [gj, and the properties for 
it/4 < 9 < 7r/2 and 9 ~ — 7r/4 are not totally agreed 
upon yet. 
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FIG. 1. Phase diagram of the general spin 1 model of Eq. 
ih after Ref. [2]. 




Hopelessly, it has not been possible until very re- 
cently to explore this phase diagram experimentally due 
to the lack of systems with a sizable biquadratic ex- 
change, which is conveniently measured by the ratio 
— j3 = J2/J1 = tan 9. It has been generally assumed 
that p is always small so that the system is always in 
the Haldane gap phase. The experiments on many spin 1 
chain compounds have indeed supported this view until a 
very recent experiment on the ID, spin 1 vanadium oxide 
LiVGc20g. In that experiment, an abrupt drop in mag- 
netic susceptibility typical of a spin-Peierls transition has 
been observed at a temperature of 22 K || . This property 
is consistent with a gapless spectrum above this temper- 
ature, and it has been argued that this behavior is most 
likely due to the presence of a significant biquadratic ex- 
change interaction || . 

Although biquadratic or higher-order spin exchange in- 
teractions have been discussed in the past and shown to 
induce anomalous magnetic properties in spin S > 1 sys- 
tems, they are usually small HQ]. It is therefore of inter- 
est to study the underlying mechanism for predominant 
biquadratic interaction. 

In this paper we propose a new microscopic mecha- 
nism leading to a significant biquadratic interaction in 
spin 1 systems. We use a microscopic model to derive 
the interaction couplings Ji and J2 in Eq. (|l|) based 
on a perturbation theory, which is compared with exact 
numerical calculations. The key in our model is to in- 
clude a third atomic orbital for an electron with slightly 
higher energy than the two singly occupied lowest orbital 
states, i.e. to consider the situation where we have a 
quasi-degeneracy of the orbitals on one site. The virtual 
electron transition via the third orbital level favors ferro- 
magnetic spin interaction, which may compensate largely 
the antiferromagnetic superexchange interaction, leading 
to the predominance of the biquadratic interaction with 
[3 > 1. We believe that this mechanism may be relevant 
to LiVGe 2 06- In that case the three orbitals come from 
the tig orbitals of the vanadium ions, and the degeneracy 
is lifted due to a small distortion of the octahedra around 
the vanadium atoms [H. More generally, we believe that 
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the exchange Hamiltonian is likely to have higher order 
contributions whenever the degeneracy between the last 
occupied and first empty orbitals is only slightly lifted. 
This situation could be realized in some other transition 
metal oxides. 

Let us consider a lattice of atoms with two outer elec- 
trons per atom. The atomic orbitals of the lowest energy 
level is two-fold degenerate, labelled by indices 1 and 2. 
There is a nearby level of orbital 3, whose energy is higher 
by an amount of A. We focus on the electron interactions 
at the same atom, and denote by U (U' < U) the direct 
Coulomb repulsion between two electrons in the same 
(different) orbitals, and by Jh > the exchange interac- 
tion. The exchange interaction favors a total spin 1 state 
for two electrons on each atom, where orbitals 1 and 2 
are both singly occupied with parallel spins, and orbital 
3 is empty. This is in accordance with Hund's rule. The 
effective spin-spin coupling arises when the electron has 
a virtual transition to the neighboring atoms. Let ti m be 
the electron hopping integral from orbital I at site i to or- 
bital m at site j with the electron spin conserved. If the 
hopping integral is small compared with U', Jh and A, 
only virtual transitions are possible, and their net effect 
is to induce an intersite spin-spin coupling. 
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FIG. 2. Typical second-order exchange processes leading 
to : a) Antiferromagnetic bilinear exchange; b) Ferromagnetic 
bilinear exchange. 

The spin Hamiltonian fly is uniquely determined by 
the total spin of two neighboring sites i and j. Let E$ the 
energy of the two-site system with total spin S, with S = 
0, 1 or 2. Then J x = \{E 2 -E 1 ), and J 2 = ^ + ^ - ^. 
Treating the hopping integrals as small parameters, the 
energy E$ can be calculated within perturbation theory. 
To illustrate the essential physics, we consider the simple 
case where t\i = t\ 2 = ti3 = 0, and t 2 2 — £33- To second 
order in the hopping integrals, and in units of t^/U', the 
spin couplings are given by || 

1 2 

4 2) = 0, (2) 

where a = t 23 /t 22 , and the a's are related to the Hund's 
couplings: 01 = (1 - iff) -1 , a 2 = (1 + 2j H )~ 1 a 3 

1 + 4j ff )-\ with j H = J H /U'. From Eqn. (§), we 
see that the biquadratic interaction vanishes in the sec- 
ond order of perturbation. This explains the smallness 

of the biquadratic coupling in most systems. The first 

(2) 

term in J{ arises from hopping to an occupied orbital 



of the neighbouring site (see Fig. 2a) and corresponds 
to superexchange interaction. It is positive and favors 
antiferromagnetic alignment of the two spins. The sec- 
ond term arises from the virtual transition via the orbital 
state 3 (see Fig. 2b). Since a\ > a 2 , the contribution to 
Ji is negative and it favors ferromagnetic alignment of 
the spins. Typical virtual transitions for the first and the 
second terms in J\ are illustrated in Fig. 2. These two ex- 
change mechanisms may compensate each other, leading 
to a very small net value of J\ . In this situation, it is nec- 
essary to extend the perturbation to include the contri- 
butions from the fourth order terms. Since A << U', Jh, 
the most important contribution arises from processes in- 
volving a virtual state of spin 1 on both sites and with 
one electron in the excited state of orbital 3. Such a pro- 
cess gives an energy correction to the second order results 
by a ratio of order of 5 = t 2 /U'A. In Fig. 3 we show a 
typical example of such processes. All the other fourth 
order contributions are smaller by a factor of A/ Jh or 
A/U' , and will be neglected in the present consideration. 

1- 

r ""\ 

— i — 

\ 

FIG. 3. Typical fourth-order process leading to a bi- 
quadratic interaction. 

Summing over these fourth order terms, we find 

4 A) =~ H+(bs-lhn (3) 

where b\ = a 2 ax, b 2 = a(a 2 + as), and 63 = (a 2 — a 3 )/2 — 
a 2 a 2 /3. The correction to J\ is simply an energy shift, 
which modifies slightly the location of the transition be- 
tween ferromagnetic and antiferromagnetic bilinear inter- 
action. The biquadratic interaction first appears in the 
fourth order in perturbation theory. While the ratio (3 is 
usually small, of the order of 5, it can be significant if the 
two contributions in the second order balance out and if 
A is relatively small. From equation (||), we also note 
that J 2 < 0, which suggests such a system is always in 
the lower half plane of the phase diagram in Fig. 1. This 
also implies (3 > if J\ > 0, and (3 < if J 1 < 0. In Fig. 
4b, we plot the value of (3 as a function of a for given 
A and Jh/U'. At a < 0.5, (3 is vanishingly small. As 
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a increases, (3 increases sharply, and the divergent point 
corresponds to the transition from antiferromagnetic to 
ferromagnetic bilinear interaction. 
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FIG. 4. Exact diagonalization results for a two-site clus- 
ter, a) Lowest energy levels in units of eV. Note the 
large gap between the three lowest levels (solid line: S=0; 
long-dashed line: S=l; short-dashed line: S=2) and the next 
one (dashed-dotted line), b) Relative biquadratic interaction 
/3 as a function of £23^22 for typical parameters (£22 = -3 eV, 
A = .3 eV, U = 6 eV, Jh = U/5). Solid line: exact result; 
dashed line: perturbation result. 

As stated above, the spin couplings are related to the 
energy levels of the two spin 1 system as a function of 
the total spin. These energies can be numerically calcu- 
lated exactly. For a two-site system with four electrons 
and three orbitals at each site, we have calculated ex- 
actly the 495 energy levels for various parameters. As 
long as A is not too small as compared to the largest 
hopping integral, the low-energy part of the spectrum 
consists of three levels with total spin S = 0, S = 1 and 
S = 2 respectively (Fig. 4a). So the low energy effective 
Hamiltonian is indeed a spin I model. From the lowest 
energies E$ for total spin S = 0, 1 and 2, we have ex- 
tracted Ji and J2- Our numerical results are in good 
agreement with the perturbation theory. In particular, 
we find a transition from antiferromagnetic to ferromag- 
netic Ji as a function of a = tii/tii-, and a sharp increase 
in (3 in the same region. This transition corresponds to 
the crossings in Fig. 4a. The exact numerical results for 
(3 as function of a for typical parameters are also plotted 
in Fig. 4b. Except for the precise value of a for which 
(3 diverges, the essential features found in the exact nu- 



merical calculations are the same as in the perturbation 
theory. Let us note that the numerical results do not 
depend qualitatively on the details of the parameters as 
long as the fundamental processes are present: There is 
always a value of £23/^22 around 1 at which (3 diverges. So 
although the simplified model presented here in order to 
have compact expressions for the fourth order perturba- 
tion results (in = £12 = ^13 = 0, £33 = £22) is not exactly 
realized in LiVGe20g, we have checked numerically that 
the appropriate extension gives the same behaviour. 

The relative value of the biquadratic interaction (3 de- 
pends sensitively on the ratio of the hopping integrals, 
but not on their magnitude. Application of a pres- 
sure on the system, if it is uniformly added, may not 
change the ratio of these hopping integrals, hence the 
biquadratic interaction sensitively. A uniaxial pressure, 
on the other hand, may change the hopping integrals in 
different proportions, and could dramatically affect the 
magnetic properties of the system. This possibility would 
be worth investigating in the case of LiVGe20g. 

Finally, let us put these results in perspective. Al- 
though this mechanism would apply in any dimension, 
there is something very special about ID. In higher 
dimension, fourth-order perturbation theory produces 
other types of four-spin interactions involving spins at 
four different sites. However, in ID these processes do 
not appear as long as the hopping to further neighbours 
can be neglected, which usually is the case. So this mech- 
anism is somehow specific to ID. This might be the rea- 
son why it was overlooked in early attempts at producing 
strong biquadratic interactions because they were done 
before Haldane's conjecture and were more concerned 
with 3D magnets. 

This mechanism is also somehow related to the prob- 
lem of orbital degeneracy. It is well known that orbital 
degeneracy leads to a Hamiltonian which is not purely 
Heisenberg, but involves a pseudo-spin for the orbital 
degrees of freedom [||. For spin 1/2 systems, for which 
most of the work has been done, the bilinear interaction 
is the only candidate - higher order terms do not appear 
because they can be rewritten in terms of the bilinear 
interaction - and the situation is a clear cut: Either the 
orbitals do not play a role, and the Hamiltonian is purely 
Heisenberg, or they do, and the Hamiltonian is a spin- 
orbital model. For spin 1 systems, we showed that, on go- 
ing from a Heisenberg to a spin-orbital model by pulling 
down an empty orbital, there is an intermediate region 
where the Hamiltonian is still a pure spin Hamiltonian, 
but with a more general interaction than just bilinear due 
to a non-trivial level crossing (Fig. 4a) . This simple idea 
should be useful in other contexts as well in looking for 
general spin Hamiltonians. In that respect, we hope that 
the present results will encourage quantum chemists to 
look more closely at exchange in situations where orbitals 
are quasi-degenerate. Quantum chemistry has proven to 
be able to reproduce very accurately on the basis of ab 
initio calculations the values of the exchange integrals 
of several systems. It should thus be possible along the 
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same lines to specify more precisely the conditions under 
which the mechanism proposed in this paper will apply. 

In conclusion, we have provided a simple but efficient 
mechanism to produce significant biquadratic interac- 
tions in spin 1 chains. It is our hope that these ideas 
will ultimately lead to an experimental investigation of 
the fascinating phase diagram of the general spin 1 model 
in one dimension. 

NOTE: After completion of this project, we became 
aware of a paper by Bhatt and Yang [R. N. Bhatt and K. 
Yang, J. Appl. Phys. 83, 7231 (1998)] who considered a 
similar problem in the context of random antifcrromag- 
netic spin chains. Their model is slightly different from 
ours (they consider the case of N degenerate orbitals with 
a simplified interaction), and they restrict themselves to 
a perturbation calculation, but our results are qualita- 
tively consistent with theirs. 

ACKNOWLEDGEMENTS: We acknowledge useful 
discussions with T. Jolicocur, K. Le Hur and O. Go- 
linelli. This work was supported in part by DOE Grant 
No. DE/FG03-98ER45687. 



[1] F. D. M. Haldane, Phys. Lett. 93A, 464 (1983); F. D. M. 

Haldane, Phys. Rev. Lett. 50, 1153 (1983). 
[2] For an overview, see e.g. J.-P. Renard, L.-P. Regnault and 

M. Verdaguer, J. Physique 49, C8-1425 (1988), and K. 

Katsumata, J. Mag. Mag. Mater. 140-144, 1595 (1995), 

and references therein. 
[3] I. Affleck, T. Kennedy, E. H. Lieb, H. Tasaki, Phys. Rev. 

Lett. 59, 799 (1987); Commun. Math. Phys. 115, 477 

(1988). 

[4] See e.g. U. Schollwock, T. Jolicoeur, T. Garel, Phys. Rev. 

B 53, 3304 (1996), and references therein. 
[5] P. Millet, F. Mila, F. C. Zhang, M. Mambrini, A. B. Van 

Oosten V. A. Pashchenko, A. Sulpice, A. Stepanov, Phys. 

Rev. Lett. 83, 4176 (1999). 
[6] E. L. Nagaev, Usp. Fiz. Nauk 136, 61 (1982) [Sov. Phys. 

Usp. 25, 31 (1982)]. 
[7] One of the exceptions is the 3-dimensional compound 

EuSe, but with a mechanism different from the one dis- 
cussed in this paper. 
[8] Here we consider the case U = U' + 2.1 h, and A << 

U',J H . 

[9] K.I. Kugel and D.I. Khomskii, Usp. Fiz. Nauk 136, 621 
(1982) [Sov. Phys. Usp. 25, 231 (1982)]. 



4 



